
""RU 2 115 137 ""CI 

(51) Int. CI 6 3 5/OQ 



miSSIAN AGENCY 
FOR PATENTS AND TRADEMARKS 

ABSTRACT OF INVENTION 



(21).(22)Application; d40l614a«». 11.05.1994 
(46) DateofpuUioatian: 10.07.1998 



(71) Appficam: 

Armizanav ^GkDtaj Egarovich, 
Ctirrvylth MP«haU KTrillov^ch. 
CherBmisin Vladirrir FtBppovich, 
ArmliDrtov Aleksi^ NIkolaovteli 

(72) Invertor; Aitnlzonov Nikolai Egorovldi, 
Chmykii MOdiaa Ktrillovich, ChBremisin VlatSmir 
RippovlGfi«Arm)2Dra>v Aloltsoj NBooloeiM 

(73) Proprietor: 

Armizonov hGteola] Egorovich, 
C3tmytm MiMiaa Kirlliavfch. 
Cherembin Vladjrror RBppovich, 
Armizor>ov Meksd Nikotaovk^ 



{54>RAN6E.PINDiNG METHOD OF LOCATION AND COMPONENTS OF VECTOR OF VELOCITY OF (EJECTS 
BY RADIO SIGNALS OF SPAC8CRAFT OF SATELLrTE RADK) NAVIGATION SYSfEMS 



(67) AtxstfBCt: 

FIELD: space radio navig^kxi and 
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object to each s^ite is found by 
rneasurement of time shilU of code sequefices 
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arnJ of components of vector of velocity and 
by nrB3Surefneiit oT received Dopplsr sliifts of 
frequency with use of systeiro tracking 

^ oaniars. O^chesinet lecaiver has one driving 

^ cfianndi and 99V0ial driven ettanndi&, it 
detennines dllfensnce t^etween ranges 

^ measured t^y driven ndoeMno channels and 
range measured by driving channel. It aiso 
datermliiee djirfarence of ratoG or dhanga or 

^ ranges between rates of change of ranges 
confuted by nneasuremanls of Ooppier shifts 

<(4 of fireqifency tiy drfven chenneSe and rale of 
change or ran^ computed by change of 

^ Deppler shffi of frequency driving 
reoetving channel. Theti deterninaUon of 
double differences of rangee end double 



differences of rates of change of ranges la 
oonduded by way of nnutual sublradlon of 
diffarQRoee of rengee and dilferenoee of 
fiatos of ehanse of ranges. EFFECT: Increased 
accuracy of datamft^ation of ooofdinatee of 
position, components of vector of velocity 
of iteienrtned obped by nevigetlon ^gnale 
of spaoecnall of satellsle radio navigation 
sysiams. wfth use of radio signals ol ground 
af>d air radiation sourc5e» and with use of 
radiations of spaoecratt of other systems 
and 9lmul8(ors. 4 2 dw9 
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The invention relates to the field of space radio 
navigation and geodesy, and may be used for 
determining positional coordinates and 
components of the velocity vector of objects. 

From Shipboard satellite navigation 
complexes, by P.S. Volosov, Y. S. Dubenko et 
al., Leningrad: Sudostroyenie Press, 1976. a 
Doppler difference range finder method is known 
for detemnining positional coordinates and 
components of the velocity vector of objects 
according to navigation radio signals of 
spacecraft in satellite radio navigation systems 
(SRNS) based on measurements of the 
differences of toppcentric distances between the 
object and two positions of one and the same 
navigation spacecraft (NSC) into sequential 
moments of time. 

The practical implementation of the known 
method is the Russian Tsikada SRNS and the 
American Transit SRNS, which are first- 
generation navigation systems. In it, integration 
of the Doppler shift of radio signals received 
over a time interval AT from an artificial earth 
navigation satellite (AENS) permits 
detenmination of the number of wavelengths 
being stacked at a variety of distances from the 
phase center of the antenna of the object's 
receiver up to two AENS positions (two positions 
of the phase center of the AENS antenna): 
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where 



ti and t2 are the transmission time of the AENS 
time mark 

Ri(ti) and R2(t2) are the distances between 
phase centers of the antennas of the object and 
AENS 

c is the speed of light 

fn is the received signal frequency 

fo is the reference signal frequency 

Fn = fn ± AfM + AfMO + Afrp + Afrp + AfAP, 

where 

fw is the frequency of the signal emitted by the 
AENS 

± Afn is the frequency instability of the emitted 
signal; 

Afno, Afrp are the unknown phase shifts caused 
by propagation of signals in the ionosphere and 
troposphere 



Afrp is the unknown phase shift caused by 
gravitational forces 

Afflp are the unknown phase shifts caused by 
other factors 
f o = f M ± Af + Af o 
where 

Afo is the known constant phase shift (frequency 
standard) 

± Af is the frequency instability of the reference 
signal. 

Allowing for what has been presented, the 
expression takes the fornn 
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From the expression it is evident that the 
integral Doppler shift is determined by two 
temis. The first iem is the measurement error 
caused by conditions for radio wave 
propagation, the Earth's gravitational field, the 
frequency instability of reference generator 
emission, and other factors. They will become 
part of the navigation equation as unknowns. 
The second term is a direct measurement of the 
change in the range inclination in wavelengths of 
the reference frequency of the object being 
detected. 

The term error of the system for tracking the 
carrier (STC), which is lacking in the navigation 
equation under consideration, also enters into 
the measurement en-or of the radio navigation 
parameter (RNP). The time function being 
tracked - the carrier frequency - has non-zero 
high order derivatives. Consequently, along with 
random (noise) errors, the actual tracking 
contour with astatism of a finite order will have 
dynamic enters caused by the presence of input- 
action derivatives of a higher order than the 
order of system astatism. Reducing random 
en-ors in an automatic frequency control system 
(AFCS) requires application of a more inertial 
feedback contour with narrowing of the low- 
frequency bandpass filter transmission band. 
With this, however the STC dynamic errors 
increase, and vice versa. 

Expressing range through coordinates of a 
rectangular geocentric coordinate system, the 
navigation equation will assume the form 
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where 

xi, yi, Zi, X2, y2, Z2 are coordinates of the phase 

center of the satellite antenna at the moments of 

time t2 and ti respectively; 

xo> yo and Zq are unknown coordinates of the 

phase center of the antenna of the object being 

determined. 

As is evident, three measurements of range 
difference in four sequential positions of the 
satellite in the orbit permit coordinates of the 
object xo, yo and to be determined. In the 
process of the measurements, we must wait 
until the range to the AENS changes by a 
sufficient amount. 

The range difference method reveals its 
distinctions at such distances (bases) between 
AENS positions in orbit, when they are 
commensurate with distances with distances 
between navigation satellites and the object 
being detected. 

According to the presented inadequate 
features of the known method, there are 

enrors caused by the STC; 

errors due to frequency instability of emission 
from the navigation satellite and the reference 
generator; 

errors caused by conditions for propagation of 
radio waves In the ionosphere and troposphere 
and other factors; 

systematic and random errors; 

low precision in determining locational 
coordinates and components of the velocity 
vector of objects when using AENS at medium- 
high and high orbits. 

Also known is a range-finding method that is 
employed as a prototype. The practical 
implementation of this method is a second- 
generation SRNS, the Russian Global Orbiting 
Navigation Satellite System (GLONASS) and the 
American Global Positioning System (GPS). The 
geometric equivalent of the final algorithm of this 
method for solving the navigation problem, 
regarding the artificial earth navigation satellites 
(AENS), is to construct a totality of positional 
surfaces, the point of intersection of which is the 
desired position of the object (On-board devices 
for satellite radio navigation, under the editorship 
of V. S. Shebshayevich; Moscow, Transport 
Press, 1988). 

To solve the navigation problem, the minimum 
necessary quantity of functional dependencies 
must be equal to the number of parameters 
under assessment. Determining the positional 
coordinates of the object is reduced to solving 
the system of equations 
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where 

Ri, , R4 are the results of measurements 

of range inclinations obtained with the aid of the 
delay tracking system (DTS): 

X, y. z are object coordinates in a geometric 
rectangular coordinate system; 

Xi, yi, Zi X4, y4, Z4 are coordinates of four 

satellites transmitted in the navigation message; 

ARt is the difference between the true 
satellite-to-object range and that measured, 
caused by a shift in the object's time scale 
relative to the time scale of the AENS; 

ARi AR4 are the measurement errors 

caused by the atmosphere, ionosphere and 
other factors. 

To determine the object's positional 
coordinates, it is necessary that simultaneously 
four satellites be found in the object's field of 
vision. As a result of solving this system of 
equations, four known quantities are 
determined: three coordinates of the object 
location (x, y, z) and a correction ARj to its time 
scale (correction to hours). 

In analogous fashion, with use of the 
measurement results aided by the satellite 
navigation system, three components are 

• a ■ 

detemnined of the velocity vector y* « and 

the correction ^ to the frequency of the 
standard for the object frequency, used to form 
the time scale: 
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where 

ti hi 

^ ^ is the range alteration rate 

(radial velocities) measured with the aid of the 
satellite navigation system; 
■ • ■ 

X' are components of the object's velocity 

vector; 

■ ■ ■ are components of the vector 
xi, yi. zi 



ABrh is the unilateral bandwidth of the UPCh 
discriminator; and Ki, K2 are the constant 
parameters dependent on the technical solution 
chosen. 

Measurement of the Doppler shift is based 
on measurement of the range increment at the 
carrier frequency, using a satellite navigation 
system. 

Assessment of the precision of range 
increment measurement is determined by an 

2L 

expression for phase dispersion of a 
scheme for tracking the canier, having the fomn 



X4. y4. Z4 

of the velocity of the four satellites; 
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T is the difference between the true velocity 
and the measured one, caused by a divergence 
of the frequencies of the frequency standards of 
the AENS and object; 

A 4 are the measurement errors 

caused by radio wave propagation conditions 
and other factors. 

Range is measured in the object's apparatus 
by measuring the time interval between the time 
marks of the code received from the satellite and 
the object's local code. 

The effectiveness of this method is 
determined mainly by the noise error of RNP 
measurement, since it is just the noise error that 
limits the compensation effect of strongly 
correlated errors. Citing from On-board devices 
for satellite radio navigation, under the editorship 
of V. S. Shebshayevich; Moscow, Transport 
Press, 1988), for assessment of the noise error, 
we employ the expression 
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where 
2 

^ ' is the dispersion of measurement noise; 
A is the length of an element of the range 
finding code; 

c/No is the relation of the signal strength to the 
spectral density of noise strength at the input of 
the receiver; 

ABcc3 is the unilateral width of the DTS band; 



where 

A is the wavelength of the carrier; 

BccH is the bandwidth of the scheme for 

tracking the carrier. 

The noise error in measurements of range 
increments at the carrier frequency is practically 
an order lower than the noise error of range 
measurements using range finding codes. 

Due to differences in, for example, GLONASS 
and GPS systems for satellite radio navigation, 
the range finding method does not permit them 
to be used in concert. 

Thus, the inadequacies of the known 
method, the prototype, are: 

errors in the system for tracking of time tag 
from the signal-to-noise ratio; 

errors in the system for tracking the carrier 
from the signal-to-noise ratio; 

errors caused by radio wave propagation 
conditions in the ionosphere and troposphere 
and other factors; 

errors caused by a time scale shift of the 
object relative to the time scale of the AENS due 
to frequency instability of satellite generators 
and the object's reference generator; 

the impossibility of making joint use of radio 
emission sources in systems having different 
assignments. 

To eliminate the ionospheric lag in known 
methods, use is made of apparatus 
compensation using dual-frequency 
measurements and of compensation employing 
corrections computed according to a priori data. 

The known prototype method is 
characterized by the following totality of actions 
on received satellite radio navigation signals: 
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reception by an N-channel receiver of dual- 
frequency radio signals N of an AENS; 

determination of the range fronn the object to 
each satellite by measuring the time shifts of 
code sequences being formed by satellite 
generators relative to code sequences being 
fomied by the object's generator; 

measurement of the range increment via 
measurement of the Increments of carrier 
phases; 

determination of the coordinates of object 

position; 

determination of the components of the 
object's velocity vector. 

The purpose of the invention is to increase 
accuracy in determining positional coordinates, 
components of the velocity vector of the object 
being detected according to SRNS spacecraft 
radio signals of ground air sources of radio 
emission, and also using the radio emission of 
satellites from other systems and their 
simulators. 

This goal is achieved in that, according to 
the proposed method, in an N-channel receiver, 
one of which is the driving one, and the others 
are the driven channels, the range difference is 
determined between ranges measured by driven 
receivers and the range measured by the driving 
receiver, and also a determination of the velocity 
differences of range changes between range 
change velocities computed as per changes in 
the Doppler shifts by driven receivers, and the 
range change velocity computed as per changes 
in the Doppler shift by the driving receiver, then 
determining the dual differences of range and 
the dual differences of range change velocity by 
means of mutual subtraction from each other of 
range differences and differences in the range 
change velocities. 

The additional distinctive features of the 
proposed method are the following. 

Via the driving and receiving devices the 
range differences between the object and two 
satellite positions are detemnined by a measured 
inten/al by measuring the phase increment of 
the carrier using automatic frequency controls of 
systems for tracking the carrier navigation radio 
signals of satellites. 

The dual range differences between the 
object and two satellite positions are determined 
by a measured interval, using measurement of 
Doppler shifts detected by receivers using 
quadrature phase detectors after multiplying 
their mean values by the measured interval. 

The driving-channel receiver receives signals 
of a simulator of satellite signals. 



The signals with Doppler frequencies are 
separated by squaring the received signals, 
subsequently returning the frequencies to the 
desired ones using frequency dividers. 

The geometric interpretation of the proposed 
method is explained using the example of a 
constellation of four GLONASS satellites and 
one GPS satellite in figure 1 . 

The GPS satellite radio signal detected by the 
receiver is the driving signal, and the channel for 
reception of the GLONASS satellite receiver is 
the driven one. Accordingly, the navigation 
signals of the GLONASS satellite, and the 
satellite receiver, are the driven ones. 

In accordance with what Is presented 
above, 



A « GPS 

Z,X Z 1 Z OP3 1 OFS z i 

3#x 9 13 GPS 1 grs a i 

4.4 4 14 GPS 1 GPS 4 1 

(1> 



where 



Ar ,&R 

* * Is the difference of the measured 



of ranges between each driven GLONASS 
satellite user and between the driving GPS 
satellite user employing range finding codes; 

2,1 3.1 4,1 are dual range 
differences. 

The geometric interpretation of determining 
coordinates and velocity vector components 
according to differences in range increments 
and the dual increment differences measured 
using increments of the carrier phases, is 
explained using the example of two satellites: a 
driving satellite and a driven GLONASS satellite, 
in figure 2. 
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The points ti, t*. iz indicate positions of the 
AENS in the orbit, which are linnits of navigation 
parameter readings (the nneasured interval). 

The differences of range increments are 
written in the following manner consequently: 



The dual range increment differences will 
take the fomn 

= Mb -A^ = iR<i >-«<«. >i- 

4,X 4.1 4Z 4j. 
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The range differences in square brackets of 
equation system (1) manifest their distinctive 
features, as was shown above at such distances 
(bases) between AENS positions in orbit when 
they are commensurate with the distance 
between a navigation satellite and the object to 
be detected. In our example, the bases are 
insignificant. To comply with this condition, 
equation system (2) is transformed into an 
identical system of equations in which this 
condition is implemented: 

3*1 32 12 3111 

4,1 4Z X Z 41 XI. 

C3> 



Thus, from the system of range differences for 
the navigation satellite orbit with identical orbit 
parameters for the constellation of 5 navigation 
satellites, the one GPS is the driving one, and 
the four GLONASS ones are driven. 

The final system of equations for dual range 
differences (1) and for dual differences in range 
increments (3), expressed through coordinates 



in a geometric rectangular coordinate system 
will take the form 
for dual range differences 
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For dual differences In range increments 
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where 

x'l, x"i. /i. y"i. z'l, z"i are coordinates 



x'4. x"4. y 4. y"4, 2*4. z% 
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of driven AENS, transmitted in navigation 
messages at the moments of time ti and t2 
respectively. 

Similarly, with the use of results from 
measurements using satellite navigation 
systems, the components of the velocity vector 
are determined: 

t> •<! <<-x> •<x,->c)*<y' -y> 'Cy -y>+ 

Sji L 1 S S £ S 

L S T S * A 
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4 T 3 1 

A# 1 1 * * * ^ 
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where 

x*i, x"i. y'l. y"i, z'l. z"i are components 



X 4, x"4. y 4, y%, Z'4, Z"4 

of the AENS velocity vector transmitted in 
navigation messages at moments of time ti and 
t2 respectively. 

It is evident in analyzing systems of 
navigation equations of dual range differences 
(4), dual differences in range increments (5) and 
velocity (6) with use of the driving and driven 



AENS radio signals and conresponding receivers 
and channels, that in the equations, the 
coordinates of the driving GPS AENS are 
compensated. Also compensated are the errors, 
caused by a divergence in the time scales and 
frequencies of GPS and GLONASS relative to 
the time scales and frequencies of the object. 

If In the navigation equations of the known 
method there are errors caused by the 
ionosphere and troposphere, then their 
differences are present in equations of the 
proposed method using dual range differences. 

To ensure a high-accuracy navigation task, 
depending on the geometric factor of 
detemiining the position in space, the satellite 
position in space is selected to be such that one 
satellite is at the zenith (ensuring high accuracy 
in determining the position along the vertical), 
which the remaining satellites are in a horizontal 
plane in directions that differ from each other by 
120 - 180^ ensuring high accuracy in 
determining the position along the horizontal, 
depending on the number of satellites used. 

Thus, the proposed method, despite, for 
example, serious differences between 
GLONASS and GPS, in methods of assigning 
ephemeredes, in the arrangement of super- 
settings of service information, in the fact that 
the systems for reporting special coordinates are 
not identical, and the difference in time scales 
formed from varying frequency and time 
standards, permits them to be used jointly, 
without making them confomi as required, i.e., 
with no organizational material processing and 
processing of the mathematical system support. 

In receiving radio navigation signals from 
the GLONASS and GPS satellites in parallel or 
sequentially, using a multiplex receiver or a 
multichannel one, and also having the GPS 
satellite be the driving one In one series of 
measurements and the GLONASS satellite as 
the driven one, and vice versa in another series, 
we can determine coordinates and velocity 
vector coordinates of the object in the GPS 
coordinate-time system or in the GLONASS 
coordinate-time system, without making them 
conform. 

Joint use of the systems ensures a certain 
universality of navigation determinations, 
reliability and trustworthy observation due to 
detemiinations being compared in differing 
systems to reveal instances where one of the 
systems has a functional violation. 

By navigational support reliability what is 
meant is the navigation system's capability at 
any moment in time to provide the object with 
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information to determine location with an 
accuracy guaranteed for the operational zone. 

By authenticity what is meant is the 
navigation system's capability to reveal 
deviations in its functioning that lead to an 
impairment in accuracy of determining 
coordinates and components of the object's 
velocity vector going beyond the limits of the 
given permissible values. 

If the system of dual-difference navigation 
equations of the proposed method, using 
measurements aided by range finding codes (1) 
Is in essence a system of range difference 
equations, then the system of navigation 
equations of dual differences in range 
increments, measured with the aid of phase 
increments of the carrier at a measured interval 
(2) is a system of equations of dual range 
differences and also permits the navigation 
problem to be solved - to detenmine position 
coordinates and components of the object's 
velocity vector. Since, as was shown above, the 
measurement accuracy of phase increment 
differences at carrier frequencies is higher by an 
order than the accuracy of time shift difference 
measurements of code sequences, then the 
accuracy by which the navigation problem is 
solved using phase increments is also higher 
than the accuracy of the solution using range 
differences. 

For purposes of further increasing the 
accuracy by which the navigation problem is 
solved using phase increments at canier 
frequencies by eliminating errors from the 
measurements caused by satellite navigation 
systems, dual differences in range increments 
are produced by separation from received 
signals with frequencies equal to Doppler 
frequency differences using quadrature phase 
detectors, to the initial outputs of which a signal 
of the driving receiver arrives, and to the 
secondary outputs of which arrive signals of 
driven receiver. Then the differences in phase 
increments are determined by multiplying mean 
values of Doppler frequency differences by the 
measured interval and detennining dual 
differences of phase increments by mutually 
subtracting them. 

The information presented corresponds to 
the implementation in apparatus, a block 
diagram of which is presented in figure 3. 
Signals with Doppler frequencies are separated 
upon receipt of phase-modulated signals with 
suppressed caniers by squaring them, then 
filtering, with subsequent return of the 



frequencies to the desired ones using frequency 
dividers. 

Signals from outputs of convolution devices 
which anive in AFCS systems of satellite 
navigation systems of receivers in figure 3, in a 
synchronism mode by delays of range finding 
codes are significantly narrow-band signals - 
restored carriers, modulated by digital data. The 
ranges for changes in carrier values are 
detemnined mainly by the Doppler shift ± 50 
kHz at frequencies of a GPS or GLONASS 
satellite), while the signal spectrum width is 
determined by the digital data spectrum (« 100 
Hz). 

AFCS signals can track signals 
corresponding to only one of the two side bands, 
and consequently, have energy losses equal to 
3 dB. Therefore, connecting the device for 
separation from received navigation signals 
equal to the Doppler frequency differences of 
the proposed method of figure 3, excluding the 
secondary side bands, does not input additional 
energy losses. 

The received and transformed satellite 
navigation radio signals arriving at quadrature 
phase detectors, already carry within 
themselves frequency shifts caused by 
generator instabilities of the spacecraft and 
object, caused by radio wave propagation 
conditions in the ionosphere and troposphere, 
shifts caused by the receiver paths, and other 
factors. Therefore, in the process of separating 
oscillations with frequencies equal to the 
Doppler frequency differences from the 
proposed method, the enumerated frequency 
deviations partially compensate for each other. 
And already with the tertiary differences their 
contribution to the accuracy of the navigation 
determinations will be insignificant. 

When used for solving the navigation 
problem of phase increments, the influences of 
phase differences on the accuracy due to the 
ionosphere and troposphere for limiting points of 
the measured interval differ little, and when 
secondary differences are formed, they are 
virtually eliminated. A special distinguishing 
feature of the proposed method is that when 
measuring differences in phase increments with 
use of oscillations equal to Doppler frequency 
differences, the signal of any radiation source 
may be used as the driving signal, be it ground- 
based, air-based, or radiation from spacecraft of 
other systems. In this instance, the main 
requirement for the receiver of the object being 
detected is the capability to receive a signal and 
to transform it in such a way that it ensures that 
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the block of quadrature phase detectors 
functions will function. There is no need to know 
the coordinates of radiation sources, their time 
systems, frequency instability and frequency 
increments due to radio wave propagation. They 
are compensated for in the process of 
navigation measurements. 

The most optimal version of device 
implementation of the proposed method is the 
version when carrier signals, modulated by 
range finding codes of the simulators are used 
as the driving signal of the object's receiver. 
Simulators permit optimization of the speed of 
frequency change specifically for each type of 
navigation system, and by that same means 
ensure they will operate in optimal fashion from 
the viewpoint of obtaining a potentially possible 
accuracy in determining position coordinates 
and components of the object's velocity vector. 

The distinguishing features of the proposed 
method are: 

Receipt by N-channel receivers of navigation 
radio signals of N satellites, one of the channels 
of which is the driving one, and the others, the 
driven ones. 

Determination of the differences in range 
increments and differences in range by 
subtracting, from the measured phase 
increments, the carrier and time shifts of code 
sequences by driven receivers of phase 
increments of carrier and time shift of code 
sequences measured by the driving device. 

iVIeasurement of dual differences of range 
increments by separating signals with 
frequencies equal to Doppler frequency 
differences, received by the driving channel and 
each driven channel of the receiver, using 
quadrature phase detectors, to the first inputs of 
which driven channel signals arrive, and to the 
second inputs of which signals of the driven 
receivers arrive, and by multiplying their mean 
values by the measured interval. 

Receipt by the driving channel of the 
receiver of radio signals of ground and air 
sources of radio emission and radio emission 
from spacecraft of other systems. 

Use of driving channels of the receiver as a 
simulator signal. 

Separation of signals with Doppler 
frequencies when receiving phase modeled 
signals with suppressed carriers by squaring 
them, with filtration and with subsequent return 
of the frequencies to the desired ones using 
frequency dividers. 

Thus, the proposed method for determining 
positional coordinates and components of the 



objects' velocity vector by radio signals of SRNS 
spacecraft is novel and exhibits significant 
differences and, when used, provides a positive 
effect consisting in increasing accuracy, 
reliability and trustworthiness of navigation 
determinations of satellite and ground radio 
navigation systems. 

What is claimed is: 

1 . Range finding method for determining the 
position and vector velocity components of 
objects by radio signals from spacecraft of 
satellite radio navigation systems, in which an N- 
channet receiver placed on the object receives 
radio navigation signals from satellites to 
determine range from the objects to each 
satellite by measuring time shifts of code 
sequences formed by satellite generators 
relative to code sequences formed by the 
object's generators, and also components of the 
velocity vector by means of measuring received 
Doppler shift frequencies using systems for 
tracking the carriers, characterized in that in the 
N-channel receiver, one of which is the driving 
one, and the other the driven channels, range 
differences are determined between ranges 
measured by the driven receivers and the range 
measured by the driving receiver, and also the 
velocity differences between range change rates 
are determined by computing according to 
measured Doppler frequency shifts by the driven 
receiver and range change rates computed 
according to measurements of the Doppler 
frequency shift by the driving receiver, then dual 
range differences and dual differences in range 
change rates are determined by means of 
mutually subtracting from each other of range 
differences and differences in range change 
rates. 

2. Method according to claim 1 , 
characterized in that range differences are 
determined by the driving and the driven 
receivers between the object and two satellite 
positions, being determined by a measured 
interval by means of measuring phase 
increments of the carriers with use of phase 
automatic frequency control of tracking systems 
of the carriers of satellite radio navigation 
signals. 

3. Method according to claim 1 , 
characterized in that the dual range differences 
between the object and two satellite positions, 
being determined by a measured interval, are 
detemnined by measuring Doppler frequency 
differences received by the receiver using 
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quadrature phase detectors, having multiplied 
their mean values by the measured interval. 

4. Method accx)rding to claims 1 - 3, 
characterized in that the receiver of the driving 
channel receives signals of a simulator of 
satellite signals. 



5. Method according to claim 3, 
characterized in that the signals with Doppler 
frequencies are separated by squaring the 
signals being received, with a subsequent return 
of the frequencies to the desired ones, using 
frequency dividers. 
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